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1. Introduction 



Recently PT— invariance of one-dimensional models in Quantum Mechanics has been 
investigated by C. Bender and collaborators [l]-[3] (see also [4]) with special emphasis on the 
spectrum of the associated Hamiltonians. Since in many cases the spectrum was found to be 
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real, PT— invariance was proposed as a generalization of standard Hermiticity. However it 
soon became clear that there are simple PT— symmetric examples, for which the spectrum 
is not real and therefore alternative criteria for reality of the spectrum were explored. 

The most systematic investigation has been performed by A.Mostafazadeh [5]- [7] (see 
also [8], [9]) elaborating on the so-called pseudo- Hermiticity: 

riHr]-^ = (1) 

with T] a Hermitian invertible operator, expressed in terms of a biorthogonal basis. This type 
of approach requires a complete solution of the spectral problem. 

For non-solvable problems, it is convenient to use the intertwining relations of SUSY 
Quantum Mechanics (SUSY QM) [10], [11] to relate pairs of Hamiltonians. One Hermitian 
and one non-Hermitian Hamiltonian may be intertwined [12], [13], [14] or, in other cases, 
non-Hermitian Hamiltonians may be intertwined [15]-[16] (Sf/SV— pseudo-Hermiticity) 
with their Hermitian conjugates. Both these constructions might lead to complex models 
with real spectra. We would like to remark that pseudo-Hermiticity differs from the 

pseudo-supersymmetry of [17]. 

While one dimensional models of this kind have been investigated in many different 
frameworks like SUSY QM [12], [13], [18], [19], PT symmetry [l]-[3], [20]-[23], [19], for two 
dimensions the advance is really at the start. To our knowledge there are only the papers 
[24] and [25] , where some complex two-dimensional potentials are studied numerically, which 
are PT— symmetric and therefore are two dimensional extensions of the x'^ + igx^ potential 
[26]. 

Within SUSY QM a class of non-trivial two-dimensional models (not allowing for separa- 
tion of variables) was constructed in [2 7] -[30]. One model of this class (generalized singular 
Morse potential) was investigated [31] in detail by two novel methods: 5'C/5'y — sepeiration 
of veiriables and two-dimensional shape inveiriance. The model is partially solvable 
(see, for example [32]) or, in alternative terminology, quasi-exactly-solvable [33], [34], this 
means that only part of eigenvalues and eigenfunctions can be found. 

In Section 2 we introduce pseudo-Hermiticity with supercharges of first and 
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second order in derivatives in one and two dimensions. In Section 3 the complexification of 
two-dimensional model of [31] is implemented in the context of — pseudo-Hermiticity 

with special attention to the appearance of levels in complex conjugated pairs. In particular, 
Subsection 3.1 contains the S't/S'y— separation of variables method, and 3.2 - the complex 
form of the two-dimensional shape invariance method. 

2. SUSY QM and ^t/^F-pseudo-Hermiticity. 

For the case of Hermitian Hamiltonians the main algebraic relation of SUSY Quantum 

Mechanics [10], [11], in all possible formulations and generalizations (for example, [35]- [39] ) 
is given by intertwining relations: 

HQ+ = Q+H- (2) 
Q-H = HQ-- Q- = (3) 

for a pair of Schrodinger operators (superpartncrs): 

H^-A + V; H^-A + V. 

These relations connect eigenfunctions with the same eigenvalues (up to zero modes of Q^): 

<pE„ = Q~4>E„; 4>En = Q^'pEr,- 

For compactness we do not introduce explicitly an index associated to possible degeneracy. 

If H and H are non-Hermitian, the two intertwining relations (2) and (3) may become 
independent, and the supercharges not necessarily Hermitian conjugate. A particular 
case, considered earlier [12], H - Hermitian, and H - not, leads to the reality of the spectrum 
of H. 
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Another possibility is to examine non-Hermitian partner Hamiltonians related by what 
we call ySC/ySy— pseudo-Hermiticity: 

Q+H = H^Q+; (4) 
HQ- = Q-Hl (5) 

Since wc deal with scalar potentials, Hermitian and complex conjugations are equivalent for 
our purposes (see also [39]). The eigenstates of //"I" with eigenvalues En will be denoted by 
{4>E*)* with eigenvalues En- They are related to those of H not only by the intertwining 
but also by direct conjugation. This has been established by using the formalism of the 
biorthogonal expansion [5] -[7]. The operator Q'^ has in general zero modes becoming non- 
invertible. 

As a consequence of Eq. (4) , one can obtain a relation which can be useful for the classi- 
fication of the spectrum: 

{En-E:^)<<PEJQUE„>=0. (6) 

First, let us notice that diagonal matrix elements in the subspace of zero modes of are 
trivially zero. So, in this case (6) does not provide any restriction on the energy En, which in 
particular can be complex having no complex conjugate partner. Clearly, a non-zero value 
of the matrix element in (6) for n = m outside the subspace, considered above, implies that 
the energy En is real, while off-diagonal non-vanishing matrix element signals that complex 
energies appear in complex conjugated pairs En = E^. 

For introducing 5'f/S'y— pseudo-Hermiticity we provide an explicit exhaustive construc- 
tion in the framework of complex SUSY QM [12] in one-dimension with first order super- 
charges [15], [16]. In this case without loss of generahty: 

H = Q-Q+ + Const; = Q+Q' + Const*; 

Q+^-d + ig{x); Q- ^+d + ig{x)^-{Q+y; d=^ 

ax 

with Const an arbitrary complex number and g{x) an arbitrary real function, leading to the 
potential: 

V{x) = —g'^{x) + ig\x) + Const. 
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In this case the equation for zero modes of Q'^ has no normahzable solution, so that 
5'C/5'y— pseudo-Hermiticity is effectively equivalent to pseudo-Hermiticity (1). Connection 
with usual PT— symmetry can be established by choosing g{x) to be of even parity. Second 
order generalizations (HSUSY) of one-dimensional supercharges [12] for S't/S'y— pseudo- 
Hermiticity were recently discussed [15], [16]. 

For two-dimensional SUSY QM models solutions of the intertwining relations (2) for 
scalar Schrodinger Hamiltonians only exist with second order supercharges [27]- [31] and 
only particular solutions have been found by suitable ansatzes. For the supercharges with 
Lorentz metrics [31] 



) + Ckdk + 49+9_ + C+9_ + C_9+ + B, 



(7) 



a solution of (2) can be provided [27]- [30] by solving the system: 



a_(c_F) = 

dlF = dtF, 



(9) 
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where x± = Xi ± X2 d± = d/dx± and C± depend only on x±, respectively: 



c+ = c,-C2 = c+{x+y, c^ = Ci + C2 = c_{x_). 



The function F, solution of (9), is represented as 



F = Fi{x+ + x_) + F2{x+ - 



The potentials V{x), V{x) and the function B{x) are expressed in terms of Fi(2xi), F2(2x2) 



and C±{x±), solutions of the system (8), (9): 



V = l{C'^ + C'_) + l{Cl + C'_) + \i^F2{x+-x.)-F,{x^ + x.)y 
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The linear character of Eq.(8) in C± allows to multiply C± by the imaginary unit keeping 
Fi_2 real. This sort of complexification renders V of Eq.(lO) complex conjugate to V. Thus 
the intertwining relations (2) in this case lead automatically to S't/S'y— pseudo-Hermiticity. 
For the two-dimensional models the existence of zero modes of can not be avoided: 
actually in the class of models studied in [31] the equation for zero modes allows separation 
of variables separation of variables). Thus these zero modes can be constructed 

from normalizable solutions of two one-dimensional equations of second order (see Section 
3 of [31]). In fact, the similarity relation, which eliminates first order derivatives from the 
supercharges, is now unitary. 

The next Section will consider the complexification of the partially solvable (quasi- 
exactly-solvable) two-dimensional model studied in [31] (generahzed singular Morse po- 
tential). The spectral problem was partially solved by two methods, one based on the 
<S't/5'y— separation of variables and the second - on the shape invariance. This model is 
a natural candidate to elucidate SUSY— pscudo-Hermiticity in two dimensions because 
it is not amenable to separation of variables. Furthermore in this model pseudo- 
Hermiticity is not equivalent to pseudo-Hermiticity due to the existence of zero modes of 
and also because the spectral problem is not exactly solvable. 

We stress that in the class of models (10) the partner Hamiltonians are not [27]-[31] 
factorizable in terms of supercharges Q^. But there are symmetry operators of fourth order 
in derivatives which can be factorized: 



3. Complex two-dimensional generalized (singular) Morse potential. 



The model is defined [31] in terms of a specific choice for C± and Fi_2 in Eqs.(lO) and 



R^Q-Q+; R^Q+Q-. 



(12) 



(11): 



C+(m) = Aiaa = (7+ (a); C_(m) = Aiaa ■ coth 



C_(a); 
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-axi 



h{xi) = 7^1(2x1) = -Afexp(-2aa;i) - 2exp( 



f2{x2) = ^F2(2x2) = +A(exp{-2ax2) - 2 exp(-Q;a;2)^ ; 



V*{x;a) = V*{x;ia) = —a^a{2a + i) smb. 



ax- 



+ 



+ A 



exp(— 2q;xi) — 2exp(— axi) + exp(— 2q;x2) — 2exp(— q;x2) 



V{x] a) = V{x] ia) = —a^a{2a — i) sinh 



ax. 



+ 



+ A 



exp(— 2Q;a:i) — 2exp(— ctxi) + exp(— 2q;x2) — 2exp(— q;x2) 



(13) 



where A is an arbitrary positive constant, and a is a real parameter. Below we will use for all 
operators and functions the "hat" notation following the definitions above. We stress that 
only for real values of the parameter a the model described above satisfies SUSY— pseudo- 
Hermiticity. 

Within this complexification the supercharges Q~^{a) are Hermitian because C± = 
C^[x±) in Eq.(7) commute with (9=p. In contrast, the supercharges Q~{a) for a G R are 
Hermitian conjugate to Q+(a), but after the complexification a ^ ia they are related by 
complex conjugation: Q~{a) — {Q~^{a))* : 



ax. 



Q^{a) = 4d+d- ± Aiaad. ± 4iaa coth(— + B{a). 



(14) 



The Hamiltonian has no definite PT— symmetry, but has a x_— refiection symmetry 
Xi ^ X2 in coordinate space (permutation symmetry). The supercharges (14) are odd. 
Therefore this model has vanishing diagonal matrix elements in (6). 

In addition, the Hamiltonian has a discrete symmetry (involution): 



V{x;a) = V{x;-a + -); 
V{x;a) = V*{x;-a). 



(15) 
(16) 
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3.1. The method of 5'C/5'y— separation of variables. 

In order to apply the method [31], one has to separate variables in Eq.(14). This can 
be achieved by the transformation, which is unitary for 6 e R : 

-iba{x+ + J coth(^^)o?a;_j = (^-^ • y^^^T^j ; (17) 

Q-{0) = U{x-,b)Q-{b)U-\x-,b) =dl-dl + ^{F,{2x,) + F.2{2x2)), (18) 

where 

2-\/A 

ii = ex.p{-axi); i = l,2. 

a 

The zero modes of can be parametrized as ^: 

fln{x;a) = U{x;aj)u),^{x;a); (19) 
Unix) = exp(-^ii^)(ei6)'''i^(-^,2s„ + l;ei)F(-n,2s„ + l;e2), (20) 

where F{—n,2sn + 1;.^) is the standard degenerate (confluent) hypergeometric function, 
reducing to a polynomial for integer n, and 

VA 1 

s„ = n-->0. 21 

a 2 

Normalizable eigenfunctions "^Eki^'id) of the Hamiltonian H{x;a) can be obtained by 
linear superposition of zero modes (19) according to [31], and their eigenvalues (after com- 
plexification a —>■ ia) read: 



Ek{a) = -Aiaa + 4a a + 2ek] = -A 



a 1, 



< 0. (22) 



Since the operator U{x; a) in (19) is unitary, the condition for normalizability of eigenfunc- 
tions '^E^{x; a) now does not depend on the parameter a and is expressed by the inequality 
(21): Sn > 0. The number of normalizable zero modes is also determined by this in- 
equality. 

^As a consequence of (18), one can derive Q+(0) = {Q-{0))* = U*{b)Q+{b){U-'^{b))*. 
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Apparently, the energies (22) have nonzero imaginary part but we remind that (6) is 
trivially satisfied by the vanishing of the matrix elements < ^Em\Q^^E„ >, since we deal 
with zero modes of Q'^. 

In order to find examples of complex conjugate energies we have to explore states outside 
the linear space of zero modes of Q^. Following the procedure of [31], we construct three 
eigenfunctions 

$(^)(f; a) = Cioix; a) ■ e^^{x; a), (23) 
with energies (-E'o(a) given in (22)): 

i«(a)=io(a)+f'n«), (24) 

and 

eW(f;a) = 1^21^'^"+'^ ¥'\a) = a\2so - l)(4ia+ 1); (25) 

e(2)(f;a) = |^2|('^'^+')(;^i + ^— 7^')(«)=4«'(^o-l)(2^a + l); (26) 

e(^) (f ; a) = - ; 7^'^ («) = «2 (4ia + 2so - l) . (27) 

Here 

_ 1 1 _ 1 1 

zi- — + —; Z2- — 

C,l C.2 Q ^2 

In contrast to the case a G R, where only ^^'^\x; a) is normalizable, all three eigenfunctions 
(25) - (27) become normalizable after a — > ia if 

So > 2. (28) 

One can argue from Eq.(l) that Q"*"^*^*-* are eigenfunctions of with the eigenvalues 
£■^'^(0). As explained after (1), this means that H{a), in addition to eigenvalues (24), has 
also the complex conjugate eigenvalues: 

^(a)(g+l>«(f; a)y = (£^® (a))*(g+l>® (f ; a))*. (29) 

The condition of normalizability of all three wave functions in (29) coincides with (28). Or- 
thogonality of all these wave functions can not be secured in general due to non-Hermiticity 
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of the Hamiltonian-'^, but a pseudo-orthogonality can be derived in agreement with the for- 
mahsm [5]- [7] of the biorthogonal expansion with particular components and 

g+(<l«) 

< > ^ i^j (30) 

I g-(<l>(^))* > = i,j = 1,2,3 (31) 

< g+d«|g-6(j) > = i^j (32) 

Equations (32) can be derived by taking into account [31] that wave functions <l>*^*\ <l>'--'^ are 
eigenfuctions of the symmetry operator (12) with different eigenvalues. Equation (30) is 
proportional to (32). Eq.(31) follows from reflection symmetry considerations ($('^ are 
even and are odd) . 

Using the symmetry property (15) one can generate a new series of levels of B.{a) — 
H(—a+^) by considering each eigenfunction (and corresponding eigenvalues) from (22), (23), 
(29) with parameter shift a — > —a + |. The partner eigenfunctions with complex conjugated 
energies can also be constructed along the same line as (29). 

3.2. The shape invariance method. 

Starting from a Schrodinger equation with potential of (13) 

H{a)4)E„{x;a) = En{a)4>EAx;a) 

with (j)En{^'j ct) - arbitrary eigenfunction, taking into account that H{a + |) = H\a), we get: 

H\a)4)E^{x; a + -) = E^{a + -)0E„(f; a + -). 

^Prom the intertwining relation: 

H{a)Q-{a)^Q-{a)H\a), 
■^Due to a;_— reflection considerations one can however easily conclude that < Q~($^'')*|$^-'^ >= 0. 
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we obtain 

H{a) (Q-{a)^EM « + ^)) = K(a + (g-(a)(^s„(f ; a + |)) 

Thereby, we are at the first step [40] of a "shape invariance chain" of wave functions and 
eigenvalues for a complex value of the parameter. Notice that the definition of potentials 
now differs from that in [31] by a constant shift Aa^a?. This leads to a vanishing of TZ{a). In 
addition, we remark that our construction will contain complex values for the parameters in 
wave functions, in operators etc only in intermediate steps, but the parameter a will always 
be kept real. 

We thus can construct an additional class of levels starting from 0^^, an eigenstate of 
the kind ii e^; Q'{^^^Y etc (sec (22), (23), (29)). These states and their complex 

conjugated are additional (particular) components of the biorthogonal basis and will fulfill 
equations similar to (30), (31), (32). 

Iterating this procedure, one generates the shape invariance chain: 



E^(f;a)^ 
associated to the energy 



Q-{a)Q-{a + ^-)Q-{a + i)...Q-{a + ^^^y-^)0£„(f; a + y) 



(33) 



In particular, for - linear combination of zero modes of (see (22)) 

<t>EM « + y ) = *^"(^; « + y) (34) 

the eigenvalues are: 

E'^{a) = -a^{4ia{sn - k) + (s„ - kf + si - 4a^). 

In order to investigate the normalizability of these eigenfunctions, it is crucial to study their 
behaviour in (^1,^2) plane following different paths. As already discussed in Section 4.3. of 
[31] for the case a G R, the relevant singularities should occur for the origin, for ^2^0 and 
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for ^1 — > ^2- By using (19), (20) and (18) for supercharges and wave functions, one can study 
the suitable critical limits in the following representation for the norm of (33): 

\\^nia)\\ ~ \\U-\a) Q-(0)C/(-^)g-(0)C/(-^)..4-(0)C/(-^) Ui2a + ik)u;J, (35) 
^ V ' 

k times 

where the x dependence has been dropped for conciseness. Normalizability can be established 
for 

,„.^_l (36) 
a Z 

which is just the normalizability condition for (34). In other words the repeated application 
of Q~ does not restrict the relevant region of normalizability. 

^Prom the same condition one can estimate the number of normalizable states generated 
by successive applications of Q~ operators in (33). It depends only on the value of the 
integer part: N = — |] and does not depend on a. For example, in the case of non- 
integer — I the total number of states can be estimated to be N(N + 1)/2, in other words 
typically y4/(2a^) for large values of N. 

As a final remark, we would like to mention that analogous results for shape invariance 
chains and their normalizability can be obtained for the model [31] (before complexification) , 
though they were not explicitly discussed. In that case are interrelated by Hermitian 
conjugation, and the calculation of the norm of the chain can be performed by an explicit 
introduction of the symmetry operator R = Q~^Q~ , provided the arguments match. After 
this it is clear that results equivalent to (36) hold. 

Our main results can be summarized as follows. In the context of the notion of 
iSf/SF— pseudo-Hermiticity two methods (already studied [31] for Hermitian models) - 
separation of variables and two-dimensional shape invariance - were used to build 
explicitly a set of eigenvalues and eigenfunctions for the complex two-dimensional singular 
Morse potential. This part of the spectrum includes complex conjugated energy pairs and in 
addition non-paired complex energies for states - linear superpositions of zero modes of the 
intertwining operators. In contradistinction to one-dimensional models, pseudo-Hermiticity 
and S't/5'F-pseudo-Hermiticity are not equivalent for two-dimensional models just due to 
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the nontrivial role of the zero modes of the supercharges. For two-dimensional scalar models 
the intertwining relations can only be solved for second order supercharges, which a priori 
may have infinite number of zero modes. Therefore in the two-dimensional case the in- 
tertwined Hamiltonians are not anymore factorizable (compare with the factorizability in 
pseudo-supersymmetry of [17]), but factorizable symmetry operators (12) of fourth order in 
derivatives exist. 

In the class of models considered in [31] and here these symmetry operators cannot be 
expressed in terms of the Hamiltonians and in general signal their integrability and possible 
degeneracy of their spectra, though no direct evidence of such degeneracy was found in the 
solved part of the spectrum. We remind that in [31] the selected energy eigenfunctions which 
were explored were simultaneously eigenfunctions of the symmetry operator, so there was 
no indication for degeneracy. 
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